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We introduce a novel scheme that combines phonon-mediated quantum logic gates in trapped ions with the 
benefits of continuous dynamical decoupling. We demonstrate theoretically that a strong driving of the qubit 
decouples it from external magnetic-field noise, enhancing the fidelity of two-qubit quantum gates. Moreover, 
the scheme does not require ground-state cooling, and is inherently robust to undesired ac-Stark shifts. The 
underlying mechanism can be extended to a variety of other systems where a strong driving protects the quantum 
coherence of the qubits without compromising the two-qubit couplings. 



A quantum processor is an isolated quantum device where 
information can be stored quantum-mechanically over long 
periods of time, but also manipulated and retrieved. This 
forbids its perfect isolation, making such a device sensitive 
to the noise introduced by either external sources, or ex- 
perimental imperfections. Additionally, the interactions be- 
tween distant quantum bits (qubits), as required to perform 
quantum logic operations, are frequently achieved by aux- 
iliary (quasi)particles whose fluctuations introduce an addi- 
tional source of noise. As emphasized recently [1], one of the 
big challenges of quantum-information science is the quest for 
methods to cope with all these natural error sources, achieving 
error rates that allow fault tolerant quantum error correction. 

We address this problem in detail for a prominent architec- 
ture, namely, trapped atomic ions [2]. Among the most rele- 
vant sources of noise in this system [4], we can list the follow- 
ing: (i) thermal noise introduced by auxiliary quasiparticles 
(i.e. phonons), (/^fluctuating external magnetic fields, (Hi) 
uncompensated ac-Stark shifts due to fluctuations in the laser 
parameters, and (iv) drifts in the phases of the applied laser 
beams. There are two different strategies to overcome these 
obstacles: (a) minimize the thermal fluctuations by laser cool- 
ing [4], searching for gates operating faster than the timescale 
set by the other noise sources [5], or (b) look for schemes that 
are intrinsically robust to the different types of noise. Among 
the latter, there are certain schemes that provide partial solu- 
tions to the above noise fluctuations, such as gates that are ro- 
bust to the thermal motion of the ions [6, 7], or the encoding 
in magnetic-field insensitive states [8] and decoherence-free 
subspaces [9]. Recently, there has been a growing effort to 
implement microwave-based quantum-information process- 
ing [11, 22, 23], exploiting the excellent control over the phase 
and amplitude of microwaves as compared to laser fields. De- 
spite these efforts, it remains a key challenge to suppress all of 
the above sources of noise. Here, we propose to accomplish 
such a step, achieving fault- tolerant error bounds, by a con- 
tinuous version of dynamical decoupling at reach of current 
technology. 

While pulsed dynamical decoupling is a well-developed 
technique [13] that has already been demonstrated for 
ions [14], its optimal combination with two-qubit gates re- 
quires a considerable additional effort [15]. Hence, simpler 
protocols are a subject of recent interest [11, 16]. We hereby 
present a decoupling scheme well suited, but not limited, to 
trapped-ion experiments with three important properties: gen- 



erality, simplicity, and robustness. This scheme is sufficiently 
general to be applied to any type ion qubits. It is simple since 
it combines two standard tools available in most ion trap lab- 
oratories, namely a carrier and a red-sideband excitation. In 
particular, it relies on the strong driving of the carrier transi- 
tion, which may be realized by laser beams for optical qubits, 
or by microwaves for hyperfine and Zeeman qubits. With this 
independent driving source, we improve simultaneously the 
performance and the speed of the gate, as compared to the 
light-shift gates [17]. Besides, this driving has the potential of 
simplifying certain aspects of previous gate schemes [18], and 
is responsible for the gate robustness at different levels. On 
top of decoupling the qubits from the magnetic-field noise, it 
suppresses the errors due to the thermal ion motion, and to un- 
compensated ac-Stark shifts. When focusing on hyperfine or 
Zeeman qubits, we can further benefit from the phase and am- 
plitude stability of microwaves. Alternatively, concatenated 
drivings may overcome amplitude fluctuations [19]. 

The system - We focus on 25 Mg + to exploit the benefits 
of microwave technology [1], although the scheme is also 
valid for other ion species. Let us consider two hyperfine 
states |0) , 1 1) with energy difference ©o to form our qubit (see 
Fig. 6(a) and Table 1). The ions arrange in a string along the 
axis of a linear Paul trap characterized by the radial and axial 
frequencies co x , C0 Z (Table I). The small radial vibrations yield 
a set of collective vibrational modes of frequencies ft),,, whose 
excitations are the transverse phonons a„, a\ [18]. If the qubit- 
qubit couplings are mediated by these quasiparticles [10, 21], 
the scheme becomes less sensitive to ion-heating/thermal mo- 
tion, and easier to operate within the Lamb-Dicke regime. 

As shown in Fig. 6(a), a pair of laser beams in a Raman con- 
figuration (red arrows) induces a transition between the qubit 
states via an auxiliary excited state. By setting their frequency 
beatnote ©l close to ft)o — ft),,, such that the detuning 8„ = 
ft>L — {(Oq — ft),,) is much smaller than the radial trap frequency 
(see Table I for the bare detuning <5l = ©l — (ft>o — (Ox)), one 
obtains the red-sideband excitation. In addition, we drive the 
carrier transition (blue arrow). For our particular qubit choice, 
this driving can be performed with microwave radiation of fre- 
quency ft)(j, such that the complete Hamiltonian is 

H c + H t = £ +^ m a+a ll e- i5 "' +H.c, 

(' in 

(i) 

where we have introduced the microwave Rabi frequency D.^, 
and the sideband coupling strengths |#i„| « ^l^J scale lin- 
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Figure 1. State-dependent dipole forces: (a) Hyperfine struc- 
ture of 25 Mg+. The states |0) = \F = 3,m F = 3) and |1) = |2,2) 
from the groundstate manifold 2 Sjy2 form our qubit. Two laser 
beams £2i , drive the red sideband via an off-resonant excited state, 
and a microwave fij directly couples to the transition, (b) Spin- 
dependent a x -force acting on a single trapped ion. The phonons 
associated to the states |+^},|— x ) are displaced in phase space ac- 
cording to Di + \ (At), Di_\ (At), and form different closed paths that 
lead to the geometric phases <j)± . (c) Trotterization of the combined 
<7 X and ct v forces. The (^-displacement D| + ^(Af) shall be fol- 
lowed by the two possible rj v displacements Z)| ± \ (Af) since | + v ) oc 
(|+v) — I— i | — v ) ) - Hence, the phase-space trajectory is not generally 
closed, (d) Schematic spin-echo refocusing of the cr' -force. By ap- 
plying a 7T-pulse Xf = of (grey box) at half the rj v -displacements, 
one obtains |± y ) — > |^ v ), such that the displacements D\ ± y are re- 
versed (dotted arrows), and the trajectory is refocused yielding a 
closed path with a well-defined geometric phase robust to the thermal 
fluctuations. 



early with the laser Rabi frequency Hl and Lamb-Dicke pa- 
rameter tj, such that \^i n \ <C 8„ (Table I). Here, we use the 
spin operators = |1,-)(0;|, and we work in the interaction 
picture rotating with the phonon and qubit frequencies. 

Sideband gates and thermal noise.- We first introduce an 
intuitive picture to understand the effects of the thermal noise 
in terms of the geometric phase gates [6, 7]. The spin-boson 
Hamiltonian (1) in the absence of the driving fid = can be 
expressed as the combination of two non-commuting spin- 
dependent forces. Each force aims at displacing the normal 
modes along a closed trajectory in phase space, which is de- 
termined by the particular eigenstates |±_ r ), |± v ) of the Pauli 
matrices o x ,o y [18]. In Fig. 6(b), we describe the action of 
the <X v -force on a single trapped ion. Depending on the spin 
state \± x ), the ion follows a different path in phase space. Af- 
ter returning to the starting point, the ion picks a geometric 
phase independent of the motional state, and only determined 
by the area enclosed by the trajectory. The spin dependence 



Table I. Specific values of trapped-ion setup 

<j\)/2n (0 x /2k (0-/27Z r\ \8y]/2tz £l L /2n £1^/2% Bq 
1.8 GHz 4 MHz 1 MHz 0.2 800 kHz 500 kHz 5.2 MHz 4 mT 



of these trajectories, together with the collective nature of the 
phonons, are the key ingredients of the two-qubit geometric 
phase gates for trapped ions in thermal motion. In contrast, 
for the much simpler single red-sideband, a x and a y forces 
are implemented, resulting in rotations around two orthogo- 
nal axis as shown in Fig. 6(c). In a Trotter decomposition, 
the consecutive concatenation of these orthogonal infinitesi- 
mal displacements spoils the closed character of the trajec- 
tory. Since the vibrations do not return to the initial state, the 
qubit and phonon states are generally entangled, and the gate 
becomes sensitive to thermal fluctuations. This intuitive pic- 
ture of the thermal noise is corroborated below by means of 
analytic and numerical arguments. 

For large detunings ^> £l]_r\, the lasers only excite vir- 
tually the vibrational modes, and the phonons can be adiabat- 
ically eliminated. In fact, it is the process where a phonon is 
virtually created by an ion, and then reabsorbed by a distant 
one, which gives rise to the effective couplings 

Hee^Jffofor, jff = -£±-& in 3?* n , (2) 

ij " °" 

which can be exploited to perform the phonon-mediated 
gates, or to implement a quantum simulation of the XY spin 
chain [7]. At certain instants of time, the unitary evolution 
corresponds to a SWAP gate [24], which performs the logic 
operation |1,0 ; ) «-» 10,1/) while leaving the remaining states 
unchanged. However, there is an additional process ignored 
in the above explanation that spoils the performance of the 
gate, namely, the phonon might be reabsorbed by the same 
ion. This leads to a residual spin-phonon coupling 

tfres = ££i(*)0f. B,(t) = Y.Binmaia,^' "-^' , (3) 

i ran 

where B,„ m = — \^in^* m (j- + According to this expres- 
sion, the ion resonance frequency fluctuates in time due to the 
collective motional dynamics. This effect can be interpreted 
as a local thermal noise in the limit of many ions, where it 
leads to dephasing. In Fig. 2(a), the critical effect of this term 
on the SWAP gate has been studied. We show the results of 
a numerical simulation of the time evolution under the full 
spin-phonon Hamiltonian, and compare it to the effective ide- 
alized description (2). As evidenced in this figure, the perfor- 
mance of the gate is severely modified by the thermal phonon 
ensemble. In fact, the swapping probabilities only approach 
the effective description (circles and squares) for ground-state 
cooled ions. As the mean phonon number is increased, the os- 
cillations get a stronger damping, and the generation of Bell 
states at half the SWAP periods (black arrows) deteriorates. 

Achieving robustness against thermal noise.- We now ad- 
dress the question of protecting the coherent spin dynamics 
from this thermal dephasing by switching D.^ ^ 0. Schemat- 
ically, this may be accomplished by refocusing the effects of 
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Figure 2. Robustness with respect to thermal noise: (a) Dy- 
namics of the swap probabilities Pw(t) (squares), Poi (t) (circles) for 
the effective gate (2), compared to the exact spin-phonon Hamilto- 
nian for a two-ion crystal with different mean phonon numbers h = 
{0,0.1,1,2,4} (solid lines). The phonon Hilbert space is truncated to 
^max = 20 excitations per mode, (b) Error e = 1 — & for the genera- 
tion of the Bell state I 1 ? - ) from | y/y) = 1 10) by the driven entangling 
gate as a function of the mean phonon number and setting n max — 14. 
As the driving power is increased, Sl^/m z £ {0,2,3.8,5.2}, the fi- 
delity approaches unity. We represent the gate error and compare it 
to the fault-tolerance (FT) threshold e t ~ 10~ 4 - 10~ 2 . 

one of the spin-dependent forces using a series of spin-echo 
pulses that invert the atomic state [Fig. 6(d)]. Note that the dy- 
namics of the dephasing (3) is characterized by the difference 
between normal-mode frequencies, | (On — (O m \ < 2(»? / (O x ) <C 
CO z , and can thus be considered as a low-frequency noise. Un- 
fortunately, this noise is still much faster than the typical gate 
times, and simple spin-echo techniques [ 1 2] would not suffice 
to get rid of the thermal errors. Instead of using complicated 
pulse sequences, we show below that the strong driving of the 
carrier transition |0) O |1) implements a continuous version 
of the refocusing of Fig. 6(d), providing a viable and simple 
mechanism for overcoming this noise. As will become clear 
later on, this driving is not only responsible for the decoupling 
from the thermal noise, but it also minimizes the undesired er- 
rors due to ac-Stark shifts and magnetic-field noise. 

A helpful account of the decoupling mechanism may be the 
following. In the dressed-state basis of the driving \± x )i — 
± \0i))/V2, the residual spin-phonon coupling becomes 

= ^B I ™|+x)«(-x|4« n e i ( n --(' B »- (U '"»' +H.c. (4) 

iinn 

For a strong driving strength Q.^, this term rotates very fast 
even for two vibrational modes that are close in frequency, 
and can be thus neglected in a rotating wave approximation. 
Note that we have assumed a vanishing phase of the driving, 
but the same argument applies for any other stable phase [18]. 

This simple argument has to be readdressed for a combina- 
tion of the carrier and red-sideband interactions (1), since the 
residual couplings are no longer described by Eq. (3). In order 
to show that a similar argument can still be applied, we have 
performed a polaron-type transformation that allows us to find 
the effective interaction and residual error terms to any desired 
order of perturbation theory [18]. We find that the dynamics 
is accurately described by the effective Hamiltonian 

H^^JffafaJ + ^^af, Jff=\Jff, (5) 

ij 



whereas the residual spin-phonon coupling is given by 
ffres = £ \ {& in a n ~ &*A) (cosh©,-^' - /sinh©,of ) (6) 

in 

such that © = Lm(^ma m - ^ r * m a nl )/28 m . By moving to the 
dressed-state basis, the residual term only involves transitions 
between the dressed eigenstates \+ x ) O |— x ), supplemented 
by the transformation on the phonons encoded in the different 
powers of ©,. Fortunately, all these transitions are inhibited 
due to the large energy gap between the dressed states set by 
£2(j. More precisely, in the strong driving regime D.^ 3> 28„ 
(see Table I), the leading order terms of the residual cou- 
pling (6) can be neglected in a rotating wave approximation. 

To check the correctness of this argument, we integrate 
numerically the complete Hamiltonian with both the side- 
band and the carrier terms (1), and take into account the 
thermal motion of the trapped ions. After the unitary evo- 
lution U(tf) = U(tf, \tf){afa^)U{\tf,Q), we calculate the fi- 
delity of producing the Bell state |*F - ) = (| 10) -i|01))/%/2, 
= max4(¥-|Tr ph {^ t (f f )p {/(f f )}| l I'-)} for different 
initial thermal states and driving strengths. The results dis- 
played in Fig. 2(b) demonstrate the promised decoupling from 
the thermal noise. We observe that the fidelity of the gate im- 
proves considerably with respect to the non-driven case when 
the amplitude lies beyond Q.^ 3> 28 n . For driving strengths 
in the 4-5 MHz range, the gate error lies within the fault- 
tolerance threshold e t ~ 10~ 2 — 10~ 4 [26] for states with mean 
phonon numbers h < 2 (see inset). However, we emphasize 
that a lower fidelity should be expected when experimental 
imperfections are taken into account. 

In Fig. 3, we describe the experimental steps required to 
create the desired Bell state. In the first step, the qubits are 
optically pumped to |0), and the initial state |i//n) = 1 10) is 
prepared by a 7C-pulse Xf obtained from a microwave reso- 
nant with the carrier transition (blue). Three comments are 
now in order. First, a magnetic field Bq m 4 mT needs to be 
applied to ensure that the Zeeman splitting between magnetic 
states is sufficiently large to avoid driving unwanted transi- 
tions. Note that the motional excitation of the ions can be ne- 
glected when driving microwave transitions owing to the van- 
ishing Lamb-Dicke factor. Second, either the ac-Stark shift 
from an off-resonant laser beam or a magnetic field gradient 
is required to effectively hide the second ion. Alternatively, 
one could use ion shuttling techniques [27]. Third, we ac- 
count for the worst possible scenario by considering (i) dif- 
ferent switching times of the lasers and microwaves, and (ii) 
imperfect timing with the microwave. By introducing global 
7T -pulses ZfZ? from the energy shift of an off -resonant strong 
microwave (yellow), we refocus the fast oscillations caused 
by the resonant microwave and correct the possible difference 
of switching times. In the second step, the two-qubit coupling 
is applied at t = to by switching on the laser beams responsible 
for the red sideband. Again, a refocusing pulse at t = tf/2 shall 
correct for the imperfect synchronization with the resonant 
microwave. Let us stress that these pulses may not be required 
in the case of accurate synchronization. In the final step, af- 
ter switching off the laser and microwaves at t = tf, the qubit 
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Figure 3. Scheme for the creation of entangled states: The two 

upper rows represent the circuit model for the two qubits, and the 
lower rows represent the switching of the laser beams (LB) and mi- 
crowaves (MW) responsible of driving the sideband and carrier tran- 
sitions. The initialization consists of optical pumping 3?, followed 
by a local 7T-pulse Xf = of . The global TT-pulses Z? 1 = o\ correct the 
possible synchronization errors. The two-qubit coupling is applied 
by switching on the laser beams responsible for the red-sideband 
(red), and the microwaves that yield the noise decoupling (blue). 
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Figure 4. Resilience to magnetic-field noise: (a) Exponential de- 
cay of the coherences (o" A (r)) for the initial states | Vb) = due 
to the magnetic dephasing. The blue circles represent the statistical 
average of the numerical integration of N = 5 • 10 3 time evolutions, 
and the solid line represents the fit to an exponential decay, (b) Error 
in the generation of the Bell state | C E ,_ ) from [ Vb) = as pre- 
dicted by the effective gate (25). The spin-phonon Hamiltonian in- 
corporates the additional magnetic-field noise, together with a strong 
microwave driving Sl^/ln = 5.2 MHz. The error is presented as a 
function of the dephasing Ti times. In the inset, we represent the time 
evolution of the fidelity for the typical noise dephasing time T 2 =5 
ms. 



state is measured by state-dependent fluorescence techniques. 
If the announced decoupling has worked correctly, this two- 
qubit gate should have generated the entangled Bell state I 1 ? - ) 
regardless of the phonon state. Let us emphasize that this gate 
is capable of producing the remaining Bell states by choos- 
ing different initial states [18] and, together with single-qubit 
rotations, becomes universal for quantum computation. 

Resilience to magnetic-field noise and ac-Stark shifts.- So 
far, we have neglected the effects of the environment. In stan- 
dard traps, the leading source of noise is due to environmental 
fluctuating magnetic fields, which limit the coherence times 
of magnetic -field sensitive states to milliseconds [4]. This is 
particularly important for multi-ion entangled states experi- 
encing super-decohere nee [21], and will also play a key role 
in our scheme considering that the two-qubit couplings lie in 
the / e ff/27r ~ 1 kHz regime. To study its consequences, we 
model the global magnetic-field noise by a fluctuating reso- 
nance frequency H n = j J^F(?)<7?. Here, F(t) is a stochastic 
Markov process [13, 14] characterized by a diffusion constant 
c and a correlation time T that determine the exponential decay 
of the coherences (ff A (f)) = e^'/ 7 ? with T 2 = 2/cT 2 [18]. By 
fixing these parameters, we can reproduce the experimentally 
observed T% « 5ms [Fig. 4(a)], and study its consequences on 
the two-qubit entangling gate. Notice that the evolution within 
the zero-magnetization subspace is not affected by this global 
noise, which can be interpreted as a decoherence-free sub- 
space. Hence, we have studied the fidelity of generating a Bell 
state that lies outside this subspace |«&~) = (|11) — i|00))/\/2- 

As shown below, the strong driving £2<j protects the qubit 
coherences from this magnetic noise without compromising 
the entangling gate, and may be considered as a continuous 
version [11, 31] of the so-called dynamical decoupling [13]. 
To single out the effects of the noise from those of the ther- 
mal motion, we have considered a ground-state cooled crys- 
tal, setting Q.i/2n = 5.2MHz to ensure that the results can 
be carried out to higher temperatures [Fig. 2(b)]. We evalu- 



ate numerically the fidelity of generating the Bell state |<t>~) 
by averaging over different samplings of the random noise 
[inset of Fig. 4(b)]. Due to the decoupling, the fidelity ap- 
proaches unity at the gate time ff = 0.7 ms. In the main 
panel of Fig. 4(b), we show that the gate error for shorter 
coherence times still lies below the fault-tolerance threshold, 
which demonstrates that the decoupling mechanism supports 
a stronger magnetic noise. Alternatively, this tells us that the 
gate tolerates smaller speeds, and thus lower Rabi frequen- 
cies of the Raman beams. This shall reduce even further the 
thermal error studied above, and the spontaneous scattering 
of photons due to the Raman beam configuration in Fig. 6(a). 
We have also calculated the fidelity of the quantum channel 

with respect to the desired quantum gate t/ e ff = e - ' 5 ^ ' [18], 
which also lies within the fault-tolerance threshold. 

As an additional advantage of our scheme, we note that it 
also minimizes the effects of uncompensated ac-Stark shifts 
that may be caused by fluctuations of the laser intensities. In 
the derivation of the red-sideband Hamiltonian (1), we have 
neglected the carrier transition by imposing a weak Rabi fre- 
quency of the laser beams. Nonetheless, the energy levels will 
be off-resonantly shifted due to an ac-Stark effect. In a real- 
istic implementation of the geometric phase gates, the shifts 
caused by off-resonant transitions to all possible states must 
be compensated by carefully selecting the laser intensities, 
frequencies, and polarizations [32]. However, fluctuations of 
these parameters will compromise this procedure introducing 
additional noise. In clear contrast, the effects of these energy 
shifts in our scheme are directly cancelled by the strong driv- 
ing in the same way that the energy shifts caused by magnetic- 
field fluctuations have been minimized. 

Let us finally comment on the effect of phase instabilities 
on the gate. In analogy to the conditional phase gate [7], but 
in contrast to M0lmer-S0rensen (MS) gates [6], our scheme 
does not depend on the slow drift of the laser phases. The 
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second-order process, whereby a phonon is virtually excited 
and then reabsorbed, is associated to the creation and sub- 
sequent annihilation of photons in the same pair of Raman 
beams, giving rise to the insensitivity to slow changes of the 
phase [see Eq. (2)]. In contrast, the MS scheme involves two 
different pairs of Raman beams, such that the cross talk leads 
to the phase sensitivity. Hence, our gate (25) only relies on the 
phase of the microwave, which is easier to stabilize as com- 
pared to the phase of the two pairs of Raman beams in the MS 
scheme. 

Conclusions and outlook- We have introduced a scheme 
that merges the notion of continuous dynamical decoupling 
with warm quantum gates in trapped ions. The decoupling, 
on top of reducing effects of external magnetic field noise, is 
also responsible for the gate robustness with respect to ther- 
mal fluctuations and ac-Stark shifts. The use of continuous 
dynamical decoupling, as opposed to pulse sequences, yields 
elegant schemes which are easier to analyze theoretically and 
realize experimentally. Hence, the direction of combining 
quantum gates and dynamical decoupling is very promising. 

Finally, we remark that our scheme can be applied to dif- 



ferent ion species, and note that it could also be combined 
with a M0lmer-S0rensen gate to protect it from magnetic field 
fluctuations and ac-Stark shifts. Moreover, we emphasize that 
these ideas may find an application beyond the field of trapped 
ions. Several quantum-information architectures also make 
use a bosonic data bus to couple distant qubits, whose coher- 
ence times are limited by the environment-induced dephasing. 
Representative examples are trapped atoms in cavities, super- 
conducting qubits coupled to transmission lines, color cen- 
ters in nanodiamonds coupled to nanomechanical resonators, 
donor states in semiconductors coupled to phonons, or quan- 
tum dots coupled to excitons. The decoupling from the noise 
induced by either the environment or by the bosonic mediator 
can be achieved along the lines here discussed. After the sub- 
mission of this work, we became aware of the interest of sim- 
ilar ideas for atoms in thermal cavities [33], and decoherence- 
free states in ion traps [34]. 
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SUPPLEMENTARY MATERIAL 



Hyperflne qubit.- We consider the isotope of magnesium 
25 Mg, which after being ionized to 25 Mg + can be confined in 
a linear Paul trap [ 1 ] . The lowest energy levels correspond to 
the valence electron lying in the s or p orbitals, which have 
a transition wavelength of X sp f=s 280 nm. The ground-state 
2 S i /2 is split due to hyperfine interactions into a couple of 
long-lived states with total angular momentum F = 2,3 and 
an energy difference of C0q/2k ps 1.8 GHz [Fig. 5]. Note 
that the transition F = 2 —> 3 is electric-dipole forbidden, and 
the decay rate F m 10~ 14 Hz is so slow that one can neglect 
spontaneous decay. Besides, by applying an additional mag- 
netic field, the magnetic sublevels are Zeeman split, and one 
can isolate two particular \F,mp) states to form our hyperfine 
qubit |0) = |3,3), |1) = |2,2). Hence, a collection of trapped 
ions can be described by the pseudo-spin Hamiltonain 



2 P 3 



F = 3 
5 MHz} t 



9.2 GHz 




(7) 



Figure 5. Hyperflne energy-level structure: (a) Hyperfine struc- 
ture of 25 Mg+. The states |0) = |F = 3,m F = 3) and 1 1) = [2,2) from 
the groundstate manifold 2 S[/2 form our qubit. Two laser beams 
Q.\,Q.2 drive the red sideband via an off-resonant excited state, and a 
microwave flj directly couples to the transition. 



where of = — |0/)(0,-|. These two-level systems can be 

coherently manipulated by either a laser in a two-photon Ra- 
man configuration Q.\,Q.2, or a direct microwave i2d [Fig. 5]. 
Note that the external magnetic field must be strong enough 
so that the drivings do not excite the population of undesired 
Zeeman sublevels. Considering the microwave and laser Rabi 
frequencies employed in this work, it suffices to set B w4 
mT, such that consecutive Zeeman sublevels are split by 20 
MHz. Alternatively, one can exploit the polarization of the 
microwave to select only the desired transition. 

Transverse phonons.- At low temperatures, the ion trap- 
ping forces balance the Coulomb repulsion, and the ions self- 
assemble in an ordered structure. The equilibrium positions 
are given by the minima of the confining and Coulomb poten- 
tials, and follow from the following system of equations 

Here, we have assumed that the transverse confinement is 
much tighter than the axial one, so that the ions arrange along 
a one-dimensional string [2]. For the two-ion crystals dis- 
cussed in the main text, we have 6 {—0.63,0.63}. Note 
that for the trapping frequencies here considered co x = 4ft), = 
27r(lMHz), the ion spacing lies in the jim range. 

We now describe the properties of the transverse collective 
modes of the ion chain around these equilibrium positions. 
The small vibrations along the jc-axis, Ax,, are coupled via the 
Coulomb interaction, which in the harmonic approximation is 




where ^ = mO) 2 (|z°.|- 3 - fyE,* l^" 3 ). and z° = 

This quadratic Hamiltonian can be diagonalized by means of 



an orthogonal transformation 




(10) 



where we have set ft = 1. Here, cv n ,a n are the phonon 
creation-annihilation operators, and \JlmGin can be in- 
terpreted as the vibrational amplitude at site r^ of the n-th 
normal mode characterized by the frequency &)„ [3]. These 
amplitudes satisfy the following relations J^m^jn = 5y> 
and Y^ij^(in%j-y^jm = where 8„ m is the Kronecker 

delta, and the normal-mode frequencies are expressed as 
co„ = cOxil+^fn) 1 ! 2 . Here, % = (co z /co x ) 2 < 1 quantifies the 
anisotropy between the axial and transverse frequencies, and 
also the width of the phonon branch co„ G [co x (l — 2<j; ) , a> v ] . 
For the two-ion crystals used in the main text co x = 4ft),, the 
diagonalization of Ijj yields <x>i = 0.968ft) r for the zig-zag 
mode i = ( 1, — 1)/ \/2, and coi — CO x for the center-off mass 
mode j^n = (1, l)/v2 Taking into account these considera- 
tions, the vibrational Hamiltonian (9), together with the qubit 
Hamiltonian, can be expressed as 

#0 = ;r£a)oCTf + £G^(4«n+2)- (U) 

L i n 

Spin-dependent dipole forces.- We consider the interac- 
tion of two laser beams in the Raman configuration shown as 
red arrows in Fig. 5. These two lasers couple the two hyperfine 
states of our qubit via an auxiliary excited state in the mani- 
fold 2 P3/2- When the transitions to the excited state F = 3 are 
far off-resonance [4], this state is seldom populated and can 
be eliminated from the dynamics. In particular, a detuning 
A/27T s» 9.2GHz exceeds both the individual Rabi frequen- 
cies Q.i,Q.2, and the decay rate from the excited state F. Ac- 
cordingly, the effective coupling between the laser beam and 



7 



the qubits is expressed as = <7+e 1 ( kL ' r ' ft)L '' + H.c. in 
the dipolar approximation. Here, we have introduced of" = 
|li)(0/|, the two-photon Rabi frequency H L = £IJ£22/2A, the 
beatnote of the laser beams <x>l = (0\ — ©2, and the effective 
laser wavevector kL = kj — k2. This wavevector is directed 
along the x-axis, and thus will only couple to the transverse 
phonon modes. In the interaction picture with respect to the 
Hamiltonian (11), the laser-ion interaction becomes 

H L = ^£ (T +e i ^^( a " e "' Bn ' +a " e, ' Bn ')e i ^-' H L)' + H.c. 

^ (12) 

By tuning the laser beatnote such that 0)l = (Oo — (O n + 8 n with 
a detuning that fulfills 8 n <C (O n , it is possible to derive the 
red-sideband Hamiltonian rj H r (t) used in the main text 
by making a Taylor expansion in the small Lamb-Dicke (LD) 
parameter r\ n = k^/^/2m(0„ <C 1, namely 

tfr(r) = £^,„(J+fl„e- i5 »' +H.c, & in = ^rj„^„. (13) 

in 

Considering the trap frequency 0) x /2n = 4 MHz, and a pair of 
orthogonal beams such that the effective wavevector is aligned 
with the x-axis, = \f2(2%/X sp ), one estimates the bare LD 
parameter to be 77 = r\ nx J (0„/ (O x ~ 0.2 (see the Table in the 
main text). Note that this derivation is only valid if the Rabi 
frequency satisfies Hl <C \(Ol — ®o| in order to neglect the 
excitation of the carrier transition in a rotating wave approx- 
imation [4] . In this work, we have considered a Raman Rabi 
frequency Q. L /27t = 500 kHz< \03l- tUo|/27T « 3.2MHz. 

Once we have derived the red-sideband term ( 1 3), let us ex- 
tend on its interpretation in terms of spin-dependent dipole 
forces outlined in the main text. By introducing the usual 
Pauli matrices = j (af ± af), the red sideband becomes 

ffr = I<%(«K + <%(0°fi d4) 
in 

where we have introduced the phonon operators 

Ol(t) = ij^e-'^'+Hc, &{ n (t) = i^,„a„e- l5 «' + H.c. 

(15) 

In order to develop an intuitive understanding of the effects 
of this Hamiltonian, let us initially consider a single ion sub- 
jected to one of the two non-commuting terms, say the a x 
part. By using a Trotter decomposition for Af — » 0, the time- 
evolution operator can be expressed as a concatenation of 
the infinitesimal unitaries U (t+At,t) « D\ + ^(At)\+ x )(+ x \ + 
Di _ x \(At)\— x}(— x\> where |± x ) stand for the eigenstates of 
a x , and D\ + j (At) are the infinitesimal displacement operators 

D {±x) (At) = e ±Aa " t= FAa*« ^ Aa = _^L T?e l5, Af . (16) 

Depending on the spin state |± A ), the ion is displaced to- 
wards a particular direction in phase space. The sequential 
application of these displacements yields a spin-dependent 
closed trajectory in phase space [Fig. 6(a)], where we set 
i^L € K- Considering the relation for the displacement op- 
erator D(a)D(J3) = e lIma ' 3 *D(a + J3), the complete time evo- 
lution will be expressed as a single displacement with a spin- 
dependent phase of a geometric origin (i.e. it only depends 




Figure 6. Spin-dependent dipole forces: (a) Spin-dependent en- 
force on a single trapped ion. The phonons associated to spin-up 
states |+.v),|— x ) are displaced in phase space according to D\ + \, 
and form different closed paths that lead to the geometric 
phases <j)±. (b) Spin-dependent a y -force, (c) Trotterization of the 
combined a x and cr v forces. The cr^-displacement D\ + \(At) shall 
be followed by the two possible a y displacements Z)| ± ^(At) since 
\+ x ) x (l+j) y})- Hence, the phase-space trajectory is not gen- 
erally closed, (d) Schematic spin-echo refocusing of the O" v -force. 
By applying a 7T-pulse Xf = af (grey box), |± y ) — > \^f y ), the dis- 
placements Di±\ are reversed (dotted arrows), such that the trajec- 
tory is refocused and the closed paths yield again a geometric phase. 

on the area enclosed by the trajectory). This is precisely what 
underlies the phase gates in ion traps [15]. By taking into ac- 
count the collective nature of the vibrational modes, one finds 
that the geometric phases depend on the spins of all ions ex- 
cited by the laser beams. This allows to implement quantum 
logic operations that do not depend on the thermal ion motion. 

If we consider the a y term alone, a similar argument would 
lead to the spin-dependent trajectories displayed in Fig. 6(b), 
and to analogous geometric phases. However, when both 
a x and a y forces are applied simultaneously, as in the red- 
sideband term (14), the Trotter decomposition leads to phase- 
space trajectories that are not closed any more [Fig. 6(c)]. 
Therefore, the geometric character of the evolution and its in- 
dependence with respect to the thermal motion of the ions are 
spoiled. Below, we show explicitly how the quantum logic 
gates will depend upon the thermal motion of the ions. 

Interactions by virtual phonon exchange.- The full 

Hamiltonian Hq + H T can be expressed in a picture where the 
phonon frequency is replaced by the detuning with respect to 
the particular sideband, namely 

H = £5„4a„ + £(^ n i f a B + ^ar4). (17) 

n in 

This expression is now amenable to perform the adiabatic 
elimination of the phonons (i.e. quasi-degenerate perturba- 
tion theory [6]) , and obtain the effective coupling between 
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|O i O i )|{n + e m }) 



+5 n 




ll^Kn-e™}) 



JVp-1 



Figure 7. Effective flip-flop interactions: In the limit of large 
detuning |J?,„| <C 5„, the energy spectrum of the Hamiltonian (17) 
is clustered in manifolds characterized by the total number of vi- 
brational excitations |{n}), such that Y.m n m = N p . The energy 
width of these manifolds is bounded by \S n — S m \ = 2t,co x , where 
£ = (co z /(D x ) 2 <C 1. The second order processes (red arrows) where 
a phonon is created and then reabsorbed (or viceversa), give rise to a 
flip-flop interaction where the qubit state is exchanged between two 
distant ions. 



the qubits. In the limit of large detuning, |J^/n| <C 8 n , only 
virtual phonon excitations can take place [Fig. 7]. There 
are two possible paths for the virtual phonon exchange be- 
tween two distant ions: i) Either the system virtually populates 
the manifold with one extra phonon, or ii) it goes through a 
lower-energy manifold with one phonon less. Since these two 
processes have an opposite detuning, their amplitudes cancel 
{&*&L)(l/8n + l/S m ) - (,^ jm )(l/8 n + \/8 m ) = 0, ex- 
cept when the exchanged phonon belongs to the same mode 
n = m. In such a case, one has to take into account the bosonic 
nature of the phonons a n a r n = 1 +a\a n , which spoils the inter- 
ference between both exchange paths, and leads to an effective 
Hamiltonian where the phonons have been eliminated 



off : 



eft' 



eff 



-1 



(18) 



At this point, we should note that there is a missing process in 
the above argument, namely that the virtually-excited phonon 
can be reabsorbed by the same ion. In this case, it does not 
lead to a phonon-mediated interaction. Due to the algebraic 
properties of the Pauli matrices, the above cancellation be- 
tween the two paths does not take place in this case, and one 
obtains an ac-Stark shift that depends on the phonons, namely 



This term has a dramatic effect on the quantum logic gates 
since it couples the spins to the phonons with a strength sim- 
ilar to that of the desired gate. Hence, the spin dynamics be- 
comes sensitive to the thermal motion of the ions, and the 
geometric character of the phase gates is lost. We note that 
this term may become a gadget, rather than a limitation, in 
the quantum simulation of Anderson localization of the vibra- 
tional excitations in an ion chain [7]. 

Polaron transformation for a strong driving.- So far, we 
have considered the phonon-mediated qubit-qubit interactions 



due to the red sideband. We now describe the effects of a 
strong driving of the carrier transition 

Hd(t) = jEjflde*^ (e 1 ^^)' + e i ( 6, » +a| «)') + H.c, (20) 

where i2j is the driving Rabi frequency, and fa its phase, and 
COd its frequency. Let us recall that for our qubit choice, this 
driving can be realized with microwaves. We consider a reso- 
nant microwave <x>i = (Oq with fi^ -C ©d, so that one can ne- 
glect the counter-rotating terms in Eq. (20). In the picture 
introduced above (17), the driven Hamiltonian becomes 

H =£ 8 n ala n + ^ £(e*a+ +H.c.) +£(^,„(T+o„ +H.c). 

n i in 

(21) 

For the sake of simplicity, we set fa = since it does not 
change the essence of the decoupling mechanism described 
below. By moving onto the dressed-state basis, |± x ); = (| 1,-) ± 
\0i))/V2, the red-sideband term becomes 



ff r(0 = I k^inil-M+xl^ - \+x)i(- 

in 

+ 1 LlM\+M+x\ - hv),KlKe- l6 »' +H.c. 



-iS„t 



(22) 



From this expression, one readily observes that the terms in- 
volving transitions between the dressed eigenstates rotate very 
fast for ^S> <5„, and their contribution to the effective inter- 
actions will be negligible. In order to give a stronger weight to 
the terms diagonal in the dressed-state basis, we perform the 
spin-dependent displacement L^*af/25„. This 

canonical transformation is formalized in terms of a Lang- 
Firsov-type polaron transformation [8] as follows 



U 



S 



L 28 n ' n 



H.c. 



(23) 



which has also been used in the context of trapped ions [10, 
23]. This unitary offers an alternative mechanism to obtain the 
spin interactions by virtual phonon exchange, and also allow 
us to calculate all the residual spin-phonon couplings to any 
desired order of the small parameter \&j n \/8 n ^ 1. In par- 
ticular, considering the algebraic properties of the spins and 
phonons, it transforms the relevant operators as follows 



UaJ t/ f = cosh(0,)af - isinh(0,)(j?, 



(24) 



where ©,■ = £„, ,^i, n a m /28 m — H.c. Therefore, it displaces the 
phonon operators, whereas it rotates the spins around the x- 
axis. These transformations yield the effective Hamiltonian 

ffeff = Y. J 'li t(J ' a j 

ij 



and the residual spin-phonon coupling 



1 reft 
4 J ij ' 



(25) 



(26) 
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To any order of perturbation theory, the residual spin-phonon 
couplings only involves terms that try to induce transitions be- 
tween the dressed eigenstates, namely oj = ;| +H.c, 
and of = |+/)(— /| +H.c Therefore, in the limit of strong 
microwave driving i2<j ^> £2l, the residual spin-phonon term 
becomes rapidly rotating and can be neglected in a rotating 
wave approximation H res w 0. For this approximation to hold, 
we have to make sure that there is no resonance with a pro- 
cess that does not conserve the number of phonons. For the 
parameters of interest to our setup, it suffices to consider that 
Cld 28„ to avoid such processes, and thus minimize the ef- 
fects of the residual spin-phonon coupling. 

By refocusing the fast rotations induced by the resonant mi- 
crowave (as considered in the main text), the effective Hamil- 
tonian becomes H e ff = £y Jfj ofo 3 -, which is responsible of 
a geometric phase gate insensitive to the thermal motion of 
the ions. After the time ff = n/(%J^f), this gate t/ e ff(ff) 
: 1 — iofo 3 -)/ \/2 generates the four entangled Bell states 



|0,0;) 

\h0j) 



72^ 
72^ 



ijOy) 



7=2^ 

■km 



(27) 



and together with single-qubit rotations it constitutes a uni- 
versal set for quantum computation. At this point, it is worth 
going back to the intuitive interpretation of the spin-dependent 
forces, and the associated geometric phase gates [Fig. 6]. The 
effect of the strong driving can be understood as a continu- 
ous version of the c x -pulses that refocuses the displacement 
caused by the o y force [Fig. 6(d)]. With this intuitive picture, 
it is clear that we have recovered the geometric character of 
the gate, and the robustness with respect to thermal motion. 

Drifts of the laser and microwave phases.- We clarify the 
role of the phases of the red-sideband and carrier excitations. 

In the red-sideband term ( 1 3), we have assumed a vanishing 
phase £2l € K. Note that the effective interaction in Eq. (18) is 
not modified by considering a different phase £2l — > Q.\ j e 1>f ' L , 
or equivalently j£; n — > ^",„e"^. Intuitively, one may con- 
sider that this phase only affects the trajectory of the displaced 
phonons, but not the geometric phase [15]. Formally, this is a 
consequence of the perturbative process where phonon excita- 
tions are virtually created and annihilated, such that the qubit- 
qubit couplings only depend on terms like &i n &*-, which are 
independent of the laser phases. According to the Raman pro- 
cess [Fig. 5], this type of terms follow from the absorption and 
emission of a photon in each of the two laser beams, thus can- 
celing the dependence on the laser phase. Note that the same 
insensitivity occurs in the presence of a strong driving (25). 
This is a clear advantage with respect to M0lmer-S0rensen 
gates [11], since those depend on the phases of two different 
sidebands that follow from a bichromatic laser beam, and are 
thus subjected to fluctuations in the optical path length. 

We have also assumed a vanishing phase of the carrier in- 
teraction (20), and argued that the essence of the decoupling 
mechanism is not altered. Let us consider the effects of a non- 



vanishing phase In this case, the dressed-state basis must 
be modified to |± d ),- = (|li)±e -i *»|0j))/V5, and thus the red- 
sideband excitation (22) in that basis becomes 



-d)i{-d\e 



-i£l d A 



-iS„t 



- d } i (-d\)a n e- ld » t +H.c. 

(28) 



Note that we can still apply the same rotating wave approx- 
imation when f2d 3> S n , and thus neglect some of the terms 
in the above expression. This encourages a modification of 
the polaron-type transformation a n — > a„ — J^i^inP °f /2<>„, 
where we have introduced the dressed spin operator of — 



e'WcT; +e~"™(T i r. By the same argument used above, we see 
that the phase appearing in the polaron transformation does 
not play any role. Hence, we would obtain an the same ef- 
fective Hamiltonian as in Eq. (25), but changing of — > of. 
We have thus arrived at the final conclusion that the effective 
Hamiltonian does only depend on the phase of the carrier ex- 
citation, but not on the phase of the red sideband. This is a 
clear advantage for our qubit choice, since the carrier phase 
comes from a microwave source, which is more stable than 
the laser beams leading to the sideband. 

Magnetic noise model.- The internal level structure of 
trapped ions can be perturbed by uncontrolled external elec- 
tric and magnetic fields. Since the Stark shifts are typically 
small [4], one is only concerned with magnetic-field fluctua- 
tions of the resonance frequency (£>q — s- COq + <9bCOo(B — Z?o) + 
5<9 B 2(Bo(B — Bo) 2 . For clock states, either there is no such 
magnetic -field dependence (e.g. |0,), |1,) have both zero mag- 
netic moment), or the linear Zeeman shift vanishes at a cer- 
tain magnetic field. In such cases, the internal state coherence 
times can reach even minutes. On the other hand, when none 
of the above conditions holds, one refers to magnetic-field 
sensitive states whose typical coherence times are reduced by 
magnetic -field fluctuations to T 2 ^l-10ms. The fluctuation 
of the resonance frequency introduces a term in the Hamil- 
tonian H n = i£,-F(f)(T?, where F(t) = —gHsB(t) fluctuates, 
jj,B is the Bohr magneton, and g the hyperfine g-factor. Low- 
frequency fluctuations (i.e. fluctuations that only take place 
between different experimental runs) can be easily prevented 
by using spin-echo sequences [12]. On the other hand, the ef- 
fects of fast-frequency fluctuations cannot be refocused unless 
complicated dynamical-decoupling sequences are used. 

We are concerned with this type of fast fluctuations, which 
lead to an exponential decay of the coherences (of(t)) as mea- 
sured by Ramsey interferometry [4]. In order to reproduce 
such an exponential decay, we use a paradigmatic model of 
noise where F(t) corresponds to a stationary, Markovian, and 
gaussian random process. Such a stochastic process is known 
as a Ornstein-Uhlenbeck process [13, 14], and is characterized 
by the following Langevin equation 

at t 

where c is the diffusion constant, T the correlation time, 
and r(f) is a gaussian white noise that fulfills (r) st = 0, 
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(r(f)r(0)) s t = 5(f). This particular stochastic differential 
equation can be integrated exactly yielding a gaussian random 
process with the following mean and variance 



(F) st =F(f )e 
(F 2 ) st -(F) 2 t = q(l- 



-('-'o)A 

-2(t-t )/x 



(30) 



which show that the correlation time T sets the time scale over 
which the process relaxes to the asymptotic values. Besides, 
the autocorrelation function (F(t)F(0)) st — c -j e ~ shows 
that T also sets the time scale such that the values of pro- 
cess are correlated or not, where we have assumed a vanishing 
mean as customary. Of primary importance to the numerical 
simulations is the update formula 



F(t 2 )=F(t l )e- 



|f(l 



(31) 



where n is a unit gaussian random variable: This formula is 
valid for an arbitrary discretization t-i = t\ + 8t, [14], and is 
thus ideally suited to for the numerical integration. 

The single-qubit coherences are affected by this type of 
noise. Considering the initial state = one obtains 



(^( f ))=±e-5( < P 2 W>-st j 



(32) 



where <p(f) ~ f df'F(f') is the random process given by the 
time integral of F(t ), and its autocorrelation function is 



((p 2 (f)) st = ct 2 ( t - x ( f- - 2e-'/ T + ie- 2 '/ T 



(33) 



In the limit of short correlation times z <C f, one obtains an 
exponential damping with a decay time T2 = 2/cx 2 . 

Let us finally comment on the modification of the Lang- 
Firsov transformation in the presence of the magnetic 
noise (31), which would lead to extra residual spin-phonon 
couplings. However, these contributions are negligible for 
the regime of interest T2 ~ 1-10 ms. This can be under- 
stood by considering the variance at long times, which leads 
to \J (F 2 ) s t = Since the magnetic-field noise in 

ion traps corresponds to the limit T <C f, we have consid- 
ered x — 0.172 throughout the text, which leads to a typ- 
ical noise strength in the 0.1-1 kHz regime. Accordingly, 
^ \J (F 2 )st, and the residual spin-phonon couplings that 
come from the Lang-Firsov transformation can be neglected. 

Fidelity of the quantum gate.- In the main text of this 
manuscript, we have quantified the efficiency of the quantum 
gate for particular entangled Bell states. However, it is also 
desirable to evaluate how close is the real time-evolution of 
the system to the desired unitary t/ e ff = e 4 ' J . This can be 
quantified by the following fidelity 



^(EW) = / dy s (y s \ul f{ g(\y s }(y s \)U eff \y s ), 



(34) 



where one integrates over the whole Hilbert space of two- 
qubit states | y/ s ), and S is the quantum channel that describes 
the real time-evolution of the trapped-ion setup 

<?(|y s )(Vs|) =Tr ph {C/(f f )(p th ® |v/ s >(v/ s |)f/(ff) t }, (35) 



where U (ff ) describes the time-evolution of the complete spin- 
phonon system including the stochastic average over the mag- 
netic field noise, and p m corresponds to the thermal Gibbs 
state for the phonons. The integral over the state space in 
Eq. (34) makes the computation of the quantum channel fi- 
delity rather inefficient, specially if one considers the addi- 
tional random sampling over the magnetic field noise. An 
alternative to this method is the evaluation of the so-called 
entanglement fidelity for the channel [24] , which can be ex- 
pressed as follows 

^e(t/eff,^) = {^\ld®Ul s S{\^ m ){^ m \)l d ®U &ii \^ m ), (36) 

where | <j) m ) — Y? a =\ \ o) <&\a) / d \s the maximally entangled 
state between the physical two-qubit system, and an ancillary 
copy of it. Hence, \a) <= {|0f,0y), |0f, ly>, \h,0j),\h, lj)}, and 
d = 4. In the expression of the entanglement fidelity, the quan- 
tum channel does not act on the ancillary qubits, namely 

&{ I M (0m I ) = Tr ph {I rf ® U (f f ) (p th <g> \<j> m ) ® U (f f ) f }. 

(37) 

The crucial point is that this entanglement fidelity is related to 
the quantum channel fidelity [25] by the following expression 



d^ e (U ei{ ,£) + l 
d+l 



(38) 



Therefore, it suffices to calculate the entanglement fi- 
delity (36), which turns out to be computationally less de- 
manding in the present setup, and then infer the value of the 
quantum channel fidelity (34) via Eq. (38). 

In Fig. 8, we present our results for this ancillary-assisted 
estimation of the quantum gate error Eaa = 1 — <^aa ( re d 
line), and compare it to a discretized version of the Haar mea- 
sure in (34), which we have denoted as £hm = 1 — #hm (yel- 
low dots). To perform the integral over state space, we sample 
randomly over A^ s = 10 3 different initial states, and average the 
corresponding fidelities. Let us also note that we have consid- 
ered A^i = 10 different dephasing times T2 € [0, 5] ms, each of 
which requires a statistical average over JV m = 5-10 3 different 
histories of the stochastic process describing the magnetic- 
field noise. Accordingly, the estimation of £hm requires the 
numerical integration of = 5 • 10 evolutions of the spin- 
phonon system, which is more demanding than the ancillary- 
based method. From the results in Fig. 8, we conclude that 
the quantum channel error approaches the fault-tolerance re- 
gion even in the presence of the magnetic-field noise. 

Comparison to state-of-the-art gate implementations.- 

Finally, let us compare our scheme for robust quantum gates 
with other existing protocols. We place a special emphasis 
on the state-of-the-art experimental implementations where 
the natural noise sources limiting the fidelities have been dis- 
cussed (see [15] and references therein). 

Let us first consider the implementations of the Cirac-Zoller 
scheme for two-ion quantum gates with the highest fidelities 
& « 0.92 for gate times ff w 0.5 ms [16]. These fidelities were 
limited by two error sources, namely, off-resonant contribu- 
tions to the carrier transition, and fluctuating laser frequencies 
and magnetic fields shifting the resonance frequency. Addi- 
tionally, these experiments require single-ion addressing and 
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Figure 8. Quantum channel error: Error e(t/ e ff) of the quantum 
channel S for the time-evolution of the two-ion system with respect 
to the desired quantum gate l/ e ff . The error is estimated according to 
the discretized version £hm of Eq. (34) (yellow dots), and compared 
to the ancillary-assisted estimate £AA( re d solid line) in (38). 



a perfect laser cooling to the motional ground-state, which 
can be hard to extend to larger ion registers. As shown in 
the main text, our scheme overcomes all these problems. The 
carrier is continuously driven so that the above off-resonant 
contribution plays no significant role (note, however, that we 
must still be careful with the off-resonant contributions to 
other Zeeman sublevels). The fluctuations of the resonance 
frequency are cured by the continuous decoupling mecha- 
nism. Besides, single-ion addressing during the gate is not 
required, and we have shown that there is no necessity of per- 
fect ground-state cooling, while potentially attaining higher 
fidelities & pa 0.99 - 0.999 with similar speeds f f pa 0.7 ms. 

Let us now compare our scheme to the so-called a z ge- 
ometric phase gates [5], which achieved a fidelity & pa 0.97 
for two-qubit gates with speeds ff pa 0.04 ms. Our gate scheme 
shares some of the advantages of this gate, such as the robust- 
ness to thermal fluctuations of the ion motion, the absence 
of single-ion addressability, and the resilience to off-resonant 
contributions to the carrier transition. It is precisely the latter 
which allows the a z -gate to attain higher speeds as compared 
to the above Cirac-Zoller gates. In fact, these speeds are re- 
quired to achieve such high fidelities in the presence of the 
errors due to due to drifts of the laser frequencies and fluctu- 
ating magnetic fields. Besides, the photon scattering due to 
the Raman beam arrangement introduces additional decoher- 
ence. Our scheme overcomes the fluctuations in the resonance 
frequency by means of the continuous decoupling, which al- 
lows us to boost the coherence times far beyond the millisec- 
ond regime, reaching higher fidelities & pa 0.99 — 0.999 for 
slower gates ff pa 0.7 ms. In fact, the dynamical decoupling 
allows still slower gates so that we may use a lower Rabi 
frequency for the Raman beams, or alternatively a larger de- 
tuning with respect to the auxiliary dipole-allowed transition, 
thus reducing further the noise by photon scattering [17]. Let 
us finally note that a source of dephasing noise common to 
both gate schemes is that of fluctuating Rabi frequencies of 
the spin-phonon couplings. This problem can be overcome by 
stabilizing the laser intensities on long time scales. 



We now compare our proposal to the M0lmer-S0rensen 
scheme for hyperfine qubits [18, 19], also known as the o§ 
geometric phase gate. The best fidelities achieved & pa 0.89 
for gate speeds ff pa 0.08 ms were limited by fluctuating ac- 
Stark shifts and photon scattering from the auxiliary excited 
state [19]. Our scheme solves the fluctuating ac-Stark shifts by 
the dynamical decoupling, and may minimize the photon scat- 
tering at the expense of using slower gates (note that the gates 
are still much faster than the coherence time enhanced via the 
decoupling mechanism). Hence, it can reach higher fidelities 
& pa 0.99 — 0.999 while preserving the advantages of the do- 
gate, which is insensitive to thermal fluctuations [18, 19], uses 
a laser beam configuration that makes it resilient to slow drifts 
of the laser phases [19], and overcomes magnetic field fluctu- 
ations by using magnetic-field insensitive states [19]. 

As emphasized in the main text, one of the challenges 
of quantum-information processing is the implementation of 
quantum logic gates below the fault-tolerance threshold. So 
far, this has only been achieved for ion optical qubits coupled 
by a M0lmer-S0rensen (MS) scheme [20], which achieved 
fidelities & pa 0.993 for two ground-state cooled ions with 
gate speeds of ff pa 0.05 ms. As a consequence of the vanish- 
ingly small spontaneous decay rate on the optical transition, 
the scheme is completely insensitive to the photon scattering 
that occurs for hyperfine qubits due to the Raman beams. The 
remaining sources of noise are caused by off-resonant contri- 
butions to the carrier transition, fluctuations of the laser fre- 
quency, and changes in the laser intensities. We note that our 
scheme can overcome two of these errors reaching similar fi- 
delities & sa 0.99 — 0.999, while it shares the sensitivity to 
fluctuations of the spin-phonon laser Rabi frequencies. The 
advantage of our decoupling protocol becomes more relevant 
as the ion number is increased, since the sensitivity of the MS 
scheme to magnetic field noise grows quadratically with the 
number of ions [21]. 

In the above schemes, most of the errors are caused by 
the particular imperfections of the laser beam arrangements. 
Recently, there has been a growing interest in developing 
quantum gates that use radio-frequency or microwave radi- 
ation [22, 23], which directly overcome the problems due 
to photon scattering in the Raman configuration, and bene- 
fit from the better control over the amplitude and phase of this 
type of radiation. In [22], two-qubit quantum gates for clock 
states with a fidelity of & pa 0.76 and speeds of ff pa 0.4 ms 
were achieved by means of oscillating magnetic-field gradi- 
ents generated in the near-field of microwave currents. In this 
case, the leading source of error corresponds to the contribu- 
tions of the oscillating gradient to the carrier transition, and 
effective ac-Zeeman shifts. In [23], a two-qubit gate with fi- 
delity J? pa 0.64 for ff pa 8 ms, has been achieved by means of 
static magnetic field gradients and radio-frequency radiation. 
We note that the leading source of noise in this case is caused 
by the decoherence induced by magnetic field fluctuations. In 
comparison, our proposal for hyperfine or Zeeman qubits con- 
stitutes an hybrid between the laser- and mircrowave-based 
methods, which can overcome some of these natural errors 
and potentially reach the fault-tolerance threshold regime. 
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